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We study response functions of integrable quantum impurity problems with an exter- 
nal field at T = using non perturbative techniques derived from the Bethe ansatz. We 
develop the first steps of the theory of excitations over the new, field dependent ground 
state, leading to renormalized (or "dressed") form- factors. We obtain exactly the low 
frequency behaviour of the dynamical susceptibility x"{ u ) i n the double well problem of 
dissipative quantum mechanics (or equivalently the anisotropic Kondo problem), and the 
low frequency behaviour of the AC noise St(u>) for tunneling between edges in fractional 
quantum Hall devices. We also obtain exactly the structure of singularities in x"{ u ) an d 
St(co). Our results differ significantly from previous perturbative approaches. 
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1. Introduction. 



One dimensional quantum impurity problems, which are of considerable physical im- 
portance, have recently been the subject of significant theoretical progresses based on their 
integrability. In particular, in a recent paper |lj , we have shown how correlation functions 
at T = and without an applied field could be obtained using massless form-factors 
|J. As a result, the response function x"(u;) in the double well problem of dissipative 
quantum mechanics || 0, and the frequency dependent conductance G(u>) for tunneling 
between edges in v = | fractional quantum Hall devices || were obtained at T = and 
without an applied field. Quantities at finite temperature and with an applied field are 
also of great interest, and potentially closer to experiments, but significantly more difficult 
to obtain, except for their DC component ||, ||. 

In this paper, we address the case T = but the applied field non-zero (this field we 
designate generically by V). We consider in particular x"( u ) i n the double well problem 
with a bias (equivalently the anisotropic Kondo problem with a field applied to the impu- 
rity), and the non-equilibrium noise St(ui) for the tunneling current in the quantum Hall 
effect in the presence of an applied voltage |L0| |llj] . 

The problem (and the origin of physically interesting phenomena) is that the applied 
field changes the structure of the ground state. The new ground state (we refer to it as 
the Fermi sea in what follows) can be obtained without much difficulty from the massless 
basis of solitons, antisolitons and breathers with factorized scattering. It is simply made 
of right moving (R) solitons and left moving (L) antisolitons filling all rapidities up to 
a Fermi value A. Excitations are obtained by adding particles and making holes in the 
sea. Due to the interacting nature of the theory, the scattering of these excitations will 
be different from the case V = 0: there is a "dressing" of the S-matrices coming from the 
Fermi sea, a phenomenon largely similar to the dressing observed in lattice regularizations 
T^j . To develop a form- factor approach for the computation of correlations, we have 



then to compute dressed form-factors: the matrix elements of physical operators between 
asymptotic states of the dressed excitations. This is a difficult problem, for technical 
reasons which will be explained later. Exact results can still be obtained however. Here, 
we study in details two questions which have been the subject of interest in the litterature: 
the low frequency behaviour of the response function in dissipative quantum mechanics and 



1 More precisely, with controlled accuracy all the way from small to large coupling. 



1 



(2) 



(1) 

( + ) 

(") □ 

1 > 

~°° P A 

Fig. 1: Multi-layered Fermi sea, only antisolitons are filled at T=0. 



11 . In addition, 



the frequency dependent shot noise in the tunnelling Hall problem |L( 
we will explore the possible existence of singularities at finite frequency in these quantities. 
The hamiltonian for the models we study here is : 



1 r° 1 

H=~ dx [87rsrn 2 + — (d z 

2 J-oo ^ng 



+ H 



B • 



This is the standard Luttinger liquid bosonic hamiltonian with a term at the boundary 
x = 0[]13[]. Both the impurity tunneling and the double well problem can be written into 
this form by standard manipulations. The boundary term depends on the problem, for 
the Hall effect it is : 



H 



B 



A 



-igVtA4>(p)/2 , AgVt-i4>(0)/2 



;i.2) 



and for dissipative quantum mechanics we have : 



H B = \ 



s #(0)/2 + g_ e -i<MP)/i 



1.3) 



In (|1.2j ) V is the voltage. The electron charge is set to e = 1, and the phases in the 
boundary (tunneling) term correspond to tunneling of Laughlin quasi-particles. In ( |1.3|) , 
e is the bias, and in this formulation, g = 1, 1/2 correspond respectively to the isotropic 
and Toulouse points. Let us recall that ( |1.3|) describes also the anisotropic Kondo problem 
with a field applied to the impurity. 

It is sometimes more appealing physically to "unfold" ( |1.2| ). Indeed, writing the field 
4> as a sum of a left and a right moving part, the integral on [— oo, 0] can be traded for 
an integral on [— oo, oo], the hamiltonian acting only on right movers. The boundary term 
then becomes an impurity term with the replacement 0(0) — > 20^(0). 
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Let us now define the quantities we will study in this language. For the impurity 
tunneling, we want to study the non-equilibrium quantum noise. In the four terminal 
geometry (see fig. 2), this noise has several components. For simplicity, we will discuss in 
this introduction the tunneling noise only. Calling It(t) the tunneling current, we consider : 



S t (u) = / dt e iut < {I t (t), I t (0)} > . (1.4) 

J — oo 

With respect to figure 2, the tunnelling current is defined as It(t) =ii?(0 + ,t) — jL(0 _ ,t). 




x=0 

Fig. 2: Four terminal system in the tunnelling v = 1/3 Hall effect. 

For dissipative quantum mechanics (or the anisotropic Kondo problem), the quantity 
we want to compute is the dynamical succeptibilityi : 



C(t) = -<[S z (t),S z (0)}> 



;i.5) 



or its Fourier transform x"( u ) ^ 0- 

As described in ]l[] the hamiltonian with the boundary interaction ( |1 . 2|) is better 
understood in the framework of the sine-Gordon theory. The idea is simple: if we add a 
term A cos <p(x) to the hamiltonian, then the resulting problem is known to be integrable, 



this is the boundary sine-Gordon model [[yj. This model is conveniently addressed in the 
basis of solitons/antisolitons, breathers, where the bulk interaction reduces to factorized 
scattering encoded in an S matrix. The boundary interaction in this basis is then described 
by a reflection matrix. In order to describe the problems we want to study, we then let 
A — ► , to get a "massless scattering" description of a free boson, which involves massless 



2 We use conventions in which h = 1 and the spins are normalised to 1. 

3 Our normalization is x"H = / ^e lLvt C(t) . 
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particles still referred to as solitons/antisolitons and breathers. Again, the boundary is 
simply described by a reflection matrix. A similar description works for ( |1.3|) . 

Having this basis, it is well known how to compute the ground state particle densities 
by using the Bethe ansatz. Knowing this ground state, a naive approach to correlations 
would be to consider matrix elements @ of the form : 



where bra and ket stand for the (shifted) ground state plus some excitations, and to com- 
pute ( |1.6|) using crossing. There are however important difficulties in this approach which 
have to do with the existence of interactions (ie a non trivial S- matrix). A more natural, 
and physically appealing, way to proceed is to think in terms of dressed or renormalised 
excitations over the ground state. 

The idea again is simple, though technically still difficult: with the help of the Bethe 
equations it is possible to compute the dressed energy and dressed scattering matrix be- 
tween excitations. Having these quantities, one can then try to write axioms for the dressed 
form- factors, in analogy with the well known case when V = 0, and then solve these 
axioms (a crucial difference is that in the present case, there is an energy scale V, and 
relativistic invariance is broken). Correlators will then follow by inserting a complete set 
of excitations between operators. 

The paper is organized as follows : In the second section, we study the new ground 
state and the structure of excitations of the auxiliary hamiltonian : 



This hamiltonian is related with the problems of interest through simple manipulations 
described in the appendix, with the correspondence e = gV for the dissipative quantum 
mechanics case. 

In the third section we discuss in considerable details the physics at low energies. We 
show that it is described by an effective free-fermion theory with renormalized parameters. 
We also discuss the physics near the thresholds of various excitation processes, laying the 
ground for further discussion of the singularities. 

Here, as in our previous works, we normalize asymptotic states such that < 6\6' >= 2tyS(9 — 



>-' 7 « < a p \0\P u ...,P q > 






e 
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The low frequency behaviour being completely under control, we are then able to 
obtain some exact results as ui — > 0. In the fourth section we discuss the behaviour of the 
dynamical succeptibility in the double well problem (or the anisotropic Kondo problem), 
and we obtain a closed expression for lim^^ x ^ as a function of V (recall e = gV, e 
the bias). We prove that Shiba's relation 



Inn = TT gx (1-8) 

u>—>0 UJ 

holds for system with bias also@. A closed expression of xo follows, which can be expanded 
as a series in y2 (i- g ) at large voltage (small coupling A), and in powers of A iX_ g at small 
voltage (large coupling). There is in particular a universal product involving small and 
large voltage properties of the static succeptibility xo '■ 

lim Xo(^)l 2(1 " 9) x [ lim V 3 - 2 °xo(V)] = F(g) (1.9) 

.V— >0 J IV— >oo 

where the function F is given explicitely in the text. 

These behaviours can be verified numerically using the framework of the Numerical 
Renormalisation Group method (NRG) [[17| which is very precise, and indeed the large field 
exponent found there is in agreement |T^| with the previous result. They could possibly also 
be checked using monte carlo simulations ||19|| , and maybe experimentally. Early results 
on x"(&), without bias, were also found in [^DJ using sum rules and scaling arguments. 

In the fifth section we discuss the noise in the tunneling problem. We obtain in 
particular the results : 

hm = — , V « T B 

lo^O UJ TV 

a ( ^ / \ \ (1-1°) 



lo^O UJ TV \V / 

Expressions for the various components of the noise in the four terminal geometry as well 
as the noise in the total current are also obtained. 

In sections 4 and 5 we also discuss the existence of potential singularities, restricting 
to 1/g an integer strictly larger than 2. We find that as well as S t (uj) should have 

singularities at all values uj = ngV, n an integer. We show that the first singularity is a 
discontinuity in the first order derivative, ie of the form \uj — gV \ (recall e = gV). We argue 
that the other singularities should be of the same form: \uj — ngV\, although we cannot 
completely prove it. All of these singularities disappear when g = 1/2 where x" is regular, 
while S t (uj) has a weaker residual singularity at u = V, of the form \u> — V| 3 . 



5 Note that this differs from the usual relation x" = ^ n 9Xo because of the normalisations 
of the spins S z to one and a factor 2tv in the definition of the Fourier transform . 
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2. The Fermi sea and the structure of excitations at T=0. 

For both problems, some simple manipulations (see the appendix) lead first to the 
consideration of the bulk hamiltonian : 

H=- dx [8ngU 2 + —(d^) 2 } + — / d t cf>, (2.1) 
2 J-oo 8tt# 4tt J_ 00 

where V is the voltage in the FQHE and e = gV for DQM (or AK). The effect of the 
impurity is then encoded into scattering matrices as explained in sections 3 and 4. 

2.1. The Fermi sea 

At T = 0, in the presence of a voltage, the ground state is factorized into a Left and 
a Right Fermi seas. Moreover, working at zero temperature allows an analytical solution 
of the Bethe equations because they become linear. Restricting to right movers, the sea 
is made of solitons filling rapidities 9 G (— oo, A] - this was discussed in J7| to which we 
refer extensively in the following. Note that all the results of this section will applly to left 
movers as well after substituting antisolitons for solitons (since the dt<p term switches sign 
in O). 

In the Fermi rapidity, A, is computed : 

e A = VG ± {0)_ 

In this formula, conventions are such that ±V/2, is the energy shift of solitons/antisolitons 
due to the voltage. We also define the "charge" of solitons and antisolitons to be ±1 i. 
The kernel G is defined by : 

G+H = M = TL?^ 1 _ e-"*. (2.3) 



where 



A = — — — rlngr+ -ln(l - a), 
2(l-s) 2 V yh 



and we define G-(u) = G+(—u). 

6 It is also convenient to define their "spin" to be ±~. Both are proportional to the conserved 
topological charge of the hamiltonian. See the appendix for more details. 
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In the following, we shall repeatedly encounter the equation : 

-A 



/(#)_ f $(e-6')f(8')de' = g(8), 

J — oo 

0e(-oo,A[; f = 0, 6>A, 



(2-4) 



with the (soliton-soliton) phase shift $(6*) = 2^^1nS , (6') : 

*(*)=/ e-™ - J 2a 7 g) U (2.5) 
J-oo 2 cosh — sinh 2( /_ g) 2tt 

As shown in 0, the function / which solves this has Fourier transform : 

where g(u>) = f_ g(9)e lu6 d9. A first example of equation ( p^l]) is provided by the density 
of solitons in the Fermi sea which satisfies (recall the convention % = 1) : 



A 

$(f)- f)')n, (Q'\dB' = . 

2tt (2.7) 



P+ (9)- f <s>(e-d') P+ (e')dd' ( 

J — oo 



0G(-oo,A[; p+ = 0, 6>A. 

Using (|2.6|) one finds : 



„ ^ _ V ^ i-) n exp [2n(l - g){9 - A - A)] 

P+U 2V5F(1-^)^ n! r(-ng)T(l-n(l-g)) ' 1 ^ 

for 6* G (— oo, A[, and = outside the sea. The density p + is discontinuous at the Fermi 
rapidity 9 = A: one has lim^^- P+(9) = t^V^Q- 

Similarly, the density of holes of solitons above the Fermi sea obeys : 

p h + (o) = [ $(o- e')p + (e')dd', e > a (2.9) 

2tt J_ 00 

This density is the analytic continuation of p+ beyond the Fermi rapidity. It is sometimes 
convenient to introduce a single, analytic quantity p{9) such that p = in the Fermi sea, 
and p = p\ outside the Fermi sea. Similar equations can be written for the densities of 
holes of antisolitons p_ and holes of breathers p^, both quantities which are defined on 
the entire rapidity interval 9 G (—00, 00). 
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2.2. Low energy excitations 

Let us now come to the low energy excitations of the theory at zero temperature. 
As we will show, when there is a voltage, there is a minimal amount of energy needed to 
create a particle or a hole: The voltage (or bias) introduces a scale in the theory, and we 
can look at low (compared with that scale) energy excitations. 

The following processes are the low energy excitations, with their associated energies 

e : 

• add a soliton: e+(6>), 9 E]A, oo) 

• destroy a soliton (or create a hole): e+(0), ^ e ( — °°5 A[; e\ < gV . 

The energy of hole and particle excitations can also be addressed analytically at zero 
temperature. For example, the excitation energy for creating a hole in the sea, e^, obeys 
the equation : 

y - e" = e\(0) - J*^ *(0 - 9 / )4(«')<W, ^ 
0G(-oo,A[; 4(#) = 0, ,9>A. 
with of course el{9) > for 9 < A and e+(A) = (this condition being actually the way 



of determining the Fermi rapidity (|2.2[) ). By differentiating with respect to 9 one sees that 
■^e\{9) = —2np + (9). Explicit solution gives : 

2 „i n! r(l-n 9 )r(§-n(l- 9 )) ^ 

the constant term in ( |2.11|) is equal to e+ (— oo), which can be obtained straightforwardly 
from ( |2.1C1| ) since in that limit this becomes a simple convolution equation. 

For 9 > A, el as defined above vanishes exactly. The excitation energy to add a 
soliton above the Fermi sea can be shown to obey : 

e+(9) = e e - \ - [ A $(9 - 9')e h + (9')d9' , (2.12) 

and thus it is the analytic continuation of — el beyond the Fermi rapidity. Using the above 
determination of el one finds explicitely : 



V V (-l) n 1 ex P 



2n(l-g) 



e + (9) = e * - -L + I-^> ] 7 ~z ~7\ " 1 ~ 77 7, (2-13) 



(A - 9 - A) 



2 4 ^ n! ™ c ^ Tr(l-g) 
n=l 



cos (n^f^j r (i - 7i/^) r (f + 



while on the other hand : 

-A 



As before, we define an analytic quantity e which is equal to — e] in the sea and e_|_ outside 
the sea. 

As an example, let us prove ( |2.12|) . Call p the density of solitons in the sea and p p 
the density of solitons above the sea (often referred to in what follows as particles). If we 
add a few solitons at rapidities 9±, . . . , 9n, we induce a change of density above the sea, 
Sp p (9p) = j- s ^ Ji 5(9 p ~ ^f)) as wen as a change of densities in the sea through the Bethe 
equations. The quantity we are looking for obeys, by definition : 

/"OO 

5E = L e+(9 p )5p p (9 p )d9 p , (2.14) 

J A 

5E = lJ A (e? - 5p(9)d9 + L J™ (e » - ^ 5p p (9 p )d9 p . (2.15) 
From the Bethe equations one has : 

/A poo 
$(9-9')5p(9')d9' + / $(0 - 9 p )5p p (9 p )d9 p . (2.16) 
-oo J A 

Let us write the solution of the general equation as : 

f(9) = I* [5(9 - 9') + L(9, 9')]g(9')d9', (2.17) 

J — oo 

where the kernel L is a symmetric function of its two arguments, and its exact expression 
is not needed in what follows. Solving ( [2.16| ) for dp and replacing in ( |2.15| ) leads to : 

e+ (9) = e 9 - j + J A $(6- 9') J [5(9' - 9") + L(9', 9")} (e e " - d9" . (2.18) 

Equation (2.12) follows then from ( ^.1U| ). 

The corresponding analysis for the momentum presents a subtlety. The reason is, that 
the momentum of excitations does not vanish at the Fermi velocity, as can easily be seen 
with the example of free fermions (g = 1/2). Therefore, one has to be very careful with 
what happens right at the Fermi surface - this difficulty did not appear for the energy 
because e vanishes at the Fermi surface anyway. We will discuss the dressed momentum 
in more details in the next section. Of course, the final result can be predicted on physical 
grounds: excitations have to be relativistic, as they obviously are in the g = 1/2 case. 
Indeed, adding a potential V amounts (in the bulk) to shifting dt<p by a constant, and this 
does not change the fact that excitations have a relativistic spectrum - more precisely, for 
excitations that have as many particles as holes, momentum = energy, while for non-neutral 
ones, momentum = energy + constant. 

In the following we define p + = e + and p\ = e+ 

While at zero temperature, there are only solitons, antisolitons as well as breathers 
do contribute to the excitation spectrum, as we now discuss. 
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2.3. Other excitations 



Other excitations occur at a finite energy above the ground state (for g < 1) and are 
obtained by the following processes : 

• add a antisoliton: e_(0), 9 6 (—00,00); e_ > (1 —g)V. 

• add an n-breather: e n (6), 9 G (—00, 00); e n > ngV 

For adding a antisoliton one finds that the excitation energy is e_(0) = V + e+(#) if 
> A and e_ (9) = V - 4(0) if < A. In particular, since the maximum energy for a hole 
is e^(— 00) = gV, we see that the threshold to add a soliton is (1 — g)V. This implies in 
particular that the low-energy processes studied above must have an energy u « (1— g)V, 
and that the limit g = 1 is highly singular in this approach. Similarly, these low energy 
processes must have an energy uj « gV, and the limit g = is also singular. 

Different physical processes can therefore occur in the excitations, which have different 
thresholds. The structure of these thresholds is quite intricate for g arbitrary. In what 
follows, we shall restrict to the usual case g = |, t an integer. Then the thresholds occur 
at energy values ?V, ie all multiples of gV. 

2.4- Charge dressing 

It is also interesting to consider the charge of excitations. Suppose we destroy a soliton 
at rapidity 9 in the sea. We define the dressed charge q+(9) = 1-j^e\ @. It is therefore 
the solution of the equation similar to e+(9) (|2.10|) , but with the left hand side y — e 9 



replaced by the opposite of the charge of the bare solitons : 

r A 

-l = q h + (9)- $(6> - 9')q h + (9')d9' '. (2.19) 
J —00 

This is actually related to other quantities we already computed in the sea: ^q+{6) = 
—2np+(9) + €+(9), from which follows the renormalized charge of the excitations at the 
Fermi surface : 

q^A) = -q+(A)- = -^/2g~, (2.20) 



where we used above values of p-(A). In (|2.20| ), q- is similarly the value of the dressed 



charge for an antisoliton added above the Fermi surface. Of course, breathers have a 
vanishing dressed charge. 



7 This definition is adequate only near the Fermi rapidity. 
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3. Physics at low energy. 

3.1. Dressed scattering 

While excitations have a relativistic dispersion relation, relativistic invariance is bro- 
ken by the existence of the energy scale e A . The next question is then, what sort of 
relativistic theory is that? To find out, we must determine the corresponding S matrix. 
First, some technicalities. Introduce the functional K such that : 

K*f(e)= ( $(e-e')f(e')de' 

J — oo 

and similarly : 

i.f(d) = f(e) 

Then ( |2.4] ) reads (/ — K) • / = g, and the general solution ( |2.17| ) reads in those terms : 

(i + L)(i-K) = i 

or equivalently (I + L)K = L. This means in turn that : 

/A pA 
[5(9 - 9') + L(9, 9')]d9' / $(6>' - 9")f(9")d9" 
-oo J —oo 

A 

L(9,9')f(9')d9'. 

-oo 

Since this is true for any function / we deduce the identity : 

/ [5(9-9')+L(9 : 9')}<S>(9' -9")d9' = L(9,9"). (3.1) 

J — oo 

Let us now consider the quantization equation allowing a few particles and holes. One 

has 

/A I* oo 

$(9-9')p(9')d9' + / $>{9-9')p p {9')d9'. (3.2) 
-oo J A 

Let us now transform this equation to make sense in the dressed theory [pi | . In this theory, 



the excitations are still particles above the sea, but they are holes in the sea. Therefore, 
the dressed equations must have on the right hand side not p, but rather ph- After some 
manipulations using the foregoing technical indentities, we obtain 

/A I* oo 

L(9,9')p h (9')d9' + 2n / L(9,9')p p (9')d9' . (3.3) 
-oo J A 

11 



The first term in this equation must be identified with the derivative of the dressed 
momentum of holes, ie one has -^p\(9) = — 27ip + (9). But it is easy to see that this 



d 



equals in turn proving as claimed in the previous section that hole excitations have 
a dispersion relation of the form momentum = energy + constant. The two other terms 
must be identified with the dressed scattering, and therefore we have i. 



' d j]nS hp (6,e') = -1(6,6'). 



(3.4) 



Similarly we have 



2in dO 
1 d 



2inde ^s PP (e,e>) = me>) 
1 d -\nS ph {e 1 e') = -L{e 1 e')., 



(3.5) 



2i7r dd' 

together with -^p + (9) = 2Tcp^_(9) = -^e + (9), proving that particle excitations are also 
relativistic. 

To determine completely the S-matrix, one has to find out the constants of integra- 
tion. This is a question that is not directly answered by the Bethe ansatz. Usually, these 
constants would be determined by using relativistic invariance, ie requiring that the full S 
matrix depends only on the difference of rapidities |12|j . Here however, relativistic invari- 
ance is broken, and the S matrix has a more complicated dependence on these rapidities. 
To compute the constants, one needs to reformulate the equivalent of crossing and unitarity 
in the dressed theory with broken relativistic invariance. Here we shall contend ourselves 
with an analysis of the low energy excitations, near the Fermi rapidity. As 9 — > A, the 
energy scale V becomes unimportant because e << V, e/V — ► 0, and relativistic invariance 
is restored. Then, L(6>, #') can be approximated by L(A,A), and integration of ([T4]) and 
( p75| ) combined with unitarity leads to the simplest solution: 

Sph(9 P , 9 h ) = -exp [2mL{A, A)(9 P - 9 h )\ 
S hp (9 h: 9 P ) = - exp [2zttL( A, A) (9 h - 9 P )] 

S PP (9 p ,9' p ) = - exp [-2mL(A,A)(9 p - 9' p )] (3.6) 
S hh (6 h , 9' h ) = - exp [-2mL(A, A)(9 h - 9' h )} 



8 In fact, the following relations hold only for 9 < 9' , since the relation between the S matrix 
and the phase shift depends on which particle has the largest rapidity 
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The overall minus sign here is necessary to ensure that no two particles or holes can 
coincide i. Let us stress that (|3.6|) is , due to the problem of integration constants, partly 
a conjecture, which we will check carefully in particular in the next subsection. For a more 
complete discussion of (p 7 !))), see the next paper 

For low energy excitations above the sea, this effective theory represented by (|3]6|) is 
simply a free fermion theory: the additional phase shifts being of the form exp[z cst(9 — 9')\ 
amount to a simple gauge transformation and can be gauged away (alternatively, form- 
factors for ( p.6| ) are free fermion form- factors, up to phases that cancel out in the correlators 
of interest). Without the sea, if one looks at say the bare S matrix of solitons, it does 
not have a well defined limit at very low energy, ie when both rapidities approach minus 
infinity, because by relativistic invariance it depends on the ratio of these energies e dl / e 02 . 
The role of the sea is to give the dressed S matrix a well defined limit at very low energy, 
ie when both rapidites approach A. It is then very natural that the limiting theory should 
be a theory of free fermions. 

In subsequent computations, it will be useful to introduce the shift function. For 
example, suppose we add a soliton above the sea at rapidity 9 P and create a hole in the 
sea at rapidity 0^. This induces a shift of the rapidities in the sea: a rapidity equal to 9 a 
initially becomes 6 a + S^ 2 W a with the conditions : 



P 

I a =^+ s(2) ^ + -L V In S(9 a - 9 P + 6®9 a - (3.7) 

P 

+ -L lnS(# a - p ) - J- In S(9 a - 6 h ). 

As usual we define the shift function by Lp+(6)S^d = F(0\6 P , 9 h ). By standard manipu- 
lations, one finds the equation obeyed by the shift: 

^-[\nS(e-e h )-\nS(e-e p )] = F(e\e p ,e h )- f $(e -e')F{e'\e p1 e h )de'. (3.8) 
^ J-00 

A formal solution of this equation follows as : 

F(6\e p ,6 h ) = [ P L(6,e')d0. (3.9) 



9 The sign of S at identical rapidities is a subject of some discussion, and largely dependent on 
the way the S matrix is defined. For us, S is the object appearing in the Fateev Zamolodchikov 
algebra. 
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3.2. The noise at low energy in a pure Luttinger liquid 

To verify the consistency of our low-energy approach, let us consider the current noise 
in a Luttinger liquid with a voltage and in the absence of impurity. Let us work within 
the effective theory of low energy excitations. In doing so, we forget completely about the 
Fermi sea and concentrate on the excitations. Let us write their energies and momentum 
e = p << V. It is convenient here to introduce new rapidities such that for particles 
e+(0 p ) = e^ p and for holes e\{9h) = e^ h . Consider then the renormalized current operator 
where the vacuum expectation value has been subtracted. For a free Fermion theory one 
has : 

< 0| : d z cp : (O,O)|/?i,02 >-+= ce^V 2 / 2 

(3.10) 

< 0| : d z <j> : (0,0)|ft, ft >+_= -ceft/V"/ 2 , 

where all the particles are assumed to be right moving, we changed notation to call a 
particle + and a hole — , and c is a constant to be determined. To fix the normalization c 
of this form-factor, let us consider the charge of the one particle state 



+ 



< 



/oo 
dx : d x (j) : {x,t)\p 2 >+ 
-oo 

/oo 
dx < 01 : d z (j) : (0, O)|0i — 2 >-+ exp\i(x — t)(e 2 — ei)] 
(3.11) 



de/d/3 



= 27^20(5(01 -0 2 ), 

where the last equality is imposed using the foregoing dressed charge computation. It 
follows that c = i\[2g. Hence at coupling g 1 all what differs from the g = 1/2 case is a 
renormalization of the two-particle form factor by y/2g. The computation of the correlator 
is then the same as for the free fermions, up to a renormalization by {^/2g) 2 since the two 
point function of the current involves the squares of form-factors. We find then the noise 
at low energy : 

S(u) * ^-\u\, (3.12) 

in agreement with the well known exact result. 

Let us stress that this little computation is actually a non trivial check. In general, 
and for vanishing voltage, the normalization of the two particle form-factor can be deter- 
mined by imposing the value of the charge. When one computes the contribution of this 
two particle form- factor to the two point function of the current, it appears that some 
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contributions are missing, due to higher form-factors: in other words higher form-factors 
are needed without voltage because the charge of excitations is 1 while the two point 
function has amplitude 2g. With a voltage, things are completely different, because the 
dressed charge is \/2g, so the (free) two-particle form-factor is sufficient to reproduce the 
two-point function amplitude. Let us stress also that there is no way one can reproduce 
the V = case by taking a "limit" V — > 0. In a massless theory, there is no such thing as 
a low energy scale. 

4. Thresholds and potential singularities 

As explained earlier, the interacting theory presents a series of thresholds : physical 
processes become allowed or forbidden when uj crosses one of the values uj = ngV. Without 
an impurity, none of these values does actually lead to a singularity in the noise, since the 
formula S = &\uj\ holds. What takes place are very special cancellations, which generally 
will be spoilt by the impurity. A simple example to see the phenomenon is the free fermions 
g = 1/2. Consider the noise close to the threshold uj = gV = -j- For uj < ^, the only 
allowed physical process is the creation of a particle hole pair. Using that the form factor 
is e^ 1 / 2 ^ 2 / 2 one can write : 



where we used that for any e\ in the interval [0, V/2], uj — e± > in the second case. Of 
course when uj > V/2, another process is possible beside exciting a particle from the sea, 
it is the creation of a pair soliton-antisoliton. The threshold for creating an antisoliton is 
at V/2, and the form factor is the same so one has then : 



and this second process adds up to the previous one to reproduce the uj dependence at all 
values of uj. Hence there is no singularity because the term involving solitons produces an 
analytic continuation of the term involving a hole. 

Since we know there is no singularity in the noise, in the absence of impurity, for any 
frequency, this means that similar cancellations must take place between various contribu- 
tions around a given threshold. For instance, (1 — g) V = ^j^V is the threshold at which 





dei = V/2, uj > V/2, 
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soliton creation is possible, while it is the maximum energy one can reach by creating 
{t — 1) holes in the sea, etc. 

Let us now investigate in more details the potential singularity at u c = gV, assuming 
g < \. At this frequency, the process involving a pair particle hole saturates while the 
process involving creation of a (first) breather above the Fermi sea kicks in. Let us consider 
in more details the process involving a breather above the Fermi sea. Once again, we will 
restrict to low energies that is to frequencies such that \u> — cu c | << u c . In that limit, the 
breathers are created at rapidities close to — oo, and in all the computations determining 
the phase shifts and the dressed S-matrix, one can replace the integrals on [—00, A] by 
integrals over the whole real axis. The renormalized breather-breather S-matrix is then 
found to be : 

! d \nS bb = + (4.1) 



2md6 °° 11 ' 1-$' 
(where the right hand side is understood in terms of Fourier transforms, $n and $1 are 
defined like $ as the logarithmic derivatives of the breather-breather and breather-soliton 
S-matrices), and this simply reproduces $n up to a rapidity renormalization. Hence, the 
renormalized scattering theory for breathers close to —00 is insensitive to the voltage, and 
behaves like an ordinary sine-Gordon model. Setting : 

e b {6) = gV + e f3 , (4.2) 

the form factor of the current can depend only on e' 3 , and by dimensional analysis, must 
simply be proportional to e 13 . Hence the leading contribution to the noise without impurity 
of the breather term is : 



A h 



/oo 
e 2l3 5 (u — u c — e 13 ) dj3 = A b (oj — iv c ), u> > uj c ; 0, u < u c (4.3) 
-00 



where A b is an amplitude we did not determine here. Similarly let us consider the particle 
hole term. Calling f {61,62) the form factor of the current between the Fermi sea and a 
state with a hole at 61 < A and an added particle at 6 2 > A, the contribution to the noise 
is proportional to : 

/A />oo 
d6 1 / ^ 2 |/(^ 1 ,^ 2 )| 2 5^-e + (^ 2 )-ei(^ 1 )]. (4.4) 
-00 J A 



Performing the 62 integral gives, for u < u> 



c 



A 1 



/ 1/(01, 02)\ 2 ^ r rd$ li 6 2 = {e + )- 1 {iV-e h + {6 1 )), (4.5) 
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while for u > u c , the integral runs from — oo to A instead, with e+(— oo) = u c . Hence the 
contributions of the particle hole term below and above uj c differ by the integral from — oo 
to (e\)~ l (u>) of the same integrand as in ((4.5|). 

Now, for iv « cu c , the values of 9\ and #2 are very far apart: 6> 2 ~ A, a finite rapidity, 
while 61 approaches the value —00. At leading order, #2) | 2 /e_|_ (6*2) must then factorize 

into a function of 6\\ f(0i), and a function of A. Setting : 

e h + {6 1 )=u c -eP\ (4.6) 

we trade 9\ for f3\. By dimensional analysis, and one has f{6{) oc e^ 1 / 2 , so introducing the 
variable x = e^ 1 , the contribution to the noise without impurity of the particle hole term 
part can be rewritten, in addition to a regular term, as : 

Aph / dx = Aph (o> c - u),u < u c (4.7) 

J OJ 

Aph being another amplitude. The two processes which "cross" at uj c behave both linearly 
in u) — u c at leading order (there are no other processes near u c , and we discuss only the 
leading singularities). Since the noise has no singularity in the absence of impurity, this 
means that (|4.3|) must be the analytic continuation of (|4.7|) , that is Ab = A p h- 

For uj c = ngV, n > 2, there are in addition background processes, that is processes 
present on both sides of u c . We shall assume that such processes do not have singular 
behaviours - equivalently that the form factors of the renormalized theory are regular 
around u c . 

Then, for uj c = 2gV, we have the one second-breather process for uj > uj c , the two 
particles - two holes process for ui < uj c : the previous analysis thus generalizes to this case. 

For uj c = ngV, n > 2, more than two processes usually cross, making the analysis more 
difficult: for instance for u = 3gV, the one third-breather process together with the three 
one-breather process are possible for uj > uj c . The absence of singularity for the noise of 
the pure system therefore does not completely constrain the amplitudes as was the case 
for uj c = gV. 

We now apply the above considerations to the two problems of physical interest. 
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5. Properties of the response function in the two-state problem. 

5.1. Low frequency behaviour. 

As a first application, we consider the low energy behaviour of the response function 
in the two-state problem with a bias V. 

In the presence of a boundary interaction we have several energy scales: Tb : uj : V. 
The function x ^ can be written as times a function F (y, -y-)- We will be able to 
use the foregoing low-energy scattering theory of excitations above the Fermi sea provided 
we restrict to uj « V: that is only the dependence of F on the first variable will be 
accessible. Observe this is enough to obtain the value lim w ^ F. 

At low energy, the only physical processes involve R-solitons and L-antisolitons. The 
boundary interaction is then fully characterized by the reflection matrix - in this 
paragraph we shall omit the isotopic indices from now on. In the presence of the boundary, 
the ground state becomes : 

||0i,-..,0 n >. (5.1) 

It is a mixture of left and right particles |](eg ||0 >= \9 > R +R\9 > L ). The rapidities are 
the same as was discussed in section 2 because there is no LR scattering. 

The low energy results are obtained by making a particle- hole pair 9 p ,9h, with the 
remaining rapidities being shifted. For the current- current correlator, this gives : 

< 9 n , ■ ■ ■ , 6>i|| Jr||6>i, 9 h , ...9 n , 9 p >< 9 P , 9 n , 9^, 9i\\ Jl\\9i, 9 n >= 
R " < 9 n , • • • , 9i\Jr\9i, 9h, ...9n, 9 P > i? --- i? x ^ 2^ 

L "' L < 6> p , n , 9h, 9~i\Jl\9\, 9 n y L - L R*{9 P ) R*{9~i) J^-R(0i) 

where the hat denotes the omitted (hole) rapidity. Observe that (|5.2| ) involves the same 
matrix elements as for the excitations in the bulk, to which the analysis of section 3 applies. 
As for the JlJl and JrJr correlators, they are not affected by the boundary interaction. 
Using the fact that R* = R~ x and inserting R*(9h)R(9h) = 1, we find that we can express 
these reflection matrices in terms of a renormalised reflection matrix for the particles and 
holes : 

R*(9 p )R(9 h ) - exp | J A [F{9\9 h ) - F(0|0 p )]^lni?(0)| R*(9 p )R(9 h ). (5.3) 
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It is then convenient to introduce a renormalized reflection matrix 



11(9) = R{9) exp 



(5.4) 



Using the same line of arguments as in [Q, one can write rerwrite x" m terms of the 
correlator of the current operator. The latter can then be expressed through form-factors. 
One finds: 



- \ 9 Re / dfadfa [n*(9 2 )n(9' 2 ) - 1] ehe&6(u - - e&' 



(5.5) 



where we have defined the 9 — > correspondence through £+(#2) = e^ 2 and e+{9 2 ) = e^ 2 . 
To clarify this, let us write the R part explicitely : 



e e * + iT B _e A + ae^ + iT B 
e e 2 — %T B e A + ae^ — %Tb 

e A G+(0) 



0-2 



where : 



a 



(5.6) 



(here we used 27rp(A) = 4§|a) and similarly : 



ae 



^ - iT, 



B 



where the additional minus sign arises from the different hole/particle parametrization. In 
the foregoing equations, Tb is a renormalized coupling, related with the bare coupling A 
by 



Xk~ 9 = 2sm(Tvg)^T(g) [ 2ge~^ G+ ^ 



29. 
An 



1-9 



GAorV • (5 - 7) 

where the expressions for G+ and A are given in section 2.1, k is a cut-off. For simplicity, 
we usually use the variable Tb in the sequel. 

The integral in ( |5.5|) can then rewritten, setting e^ 2 = x : 



/ [K{x)n*{x-u)-\]dx, 
Jo 



(5.8) 



(note the dependence 011 x — u instead of u — x for the same reason as above). Obviously 
this vanishes when uj = 0. Because \R\ 2 = 1, the term linear in 10 vanishes too. At second 



order one would have : 



to 



d 

o dx 
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InTZdx 



and this will not contribute when we take the real part because R is a pure phase. The 
imaginary part of this expression, though, contributes to the real part of the static spin-spin 
succeptibility x = x' + ^x" ■ Its expression is given by : 



lim ly'M = xo 



) a 



2^g 



(5.9) 



For x"> we find that the first non trivial term goes like w 3 as expected. Collecting all 
terms one finds : 

v"C^ a 2 ( a d \ 2 

' -" 1 - " ' . (5.10) 



u>^0 UJ 



From these two previous expressions we can prove Shiba's relation in the presence of a 
bias : 

hm = n gx , (5.11) 



which is exactly the same as the one without biaM 

In particular for g = 1/2 there is no renormalization due to the sea, 1Z(9) = R(9) 



e s -iT B 



A and thus, using y = e A in that case : 

X» 2 



lim 

oj— >0 UJ 



T 2 

B 



7T 2 (V* ~^ 2 ' 



in agreement with the exact result |23l - 

Let us rewrite ( |5.10| ) in a more explicit form : 



(5.12) 



0;— >0 OJ 



a 



-2A 



4giv z 



(coshL9 b ) + I_ 00 Te F{m 



}=A cosh(0' - 9 B ) 



d9' 



(5.13) 

Let us first investigate the behaviour as V/Tb — > 0. At leading order the term in the 
bracket can be rewritten as : 



2 



e A + 



2 

Tb 



e A + 



e 6 ' ^ F(9'\9)\ 9=A 



e d L(9 1 A)d9 



4tt 



10 This relation differs from the usual lim^^o x" { u) )/ u) = ^9X% because of a different conven- 
tion for Fourier transforms and spin normalization. 
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where we used (|2.7|) , ( |2.17| ), (^9|). Hence, at very small voltage one has : 

This result is identical with the value of lim^^o x at vanishing voltage [Q. This 
proves that the limit to — > and V —>■ commute, as seems clear from the NRG method 
|T7fl . Notice that a priori this result is non trivial, the structure of excitations being very 
different at vanishing and non-vanishing voltage in a Luttinger liquid. 

To obtain more insight in the behaviour of lim^^o x — , some transformations are 
useful, which are related with the standard Bethe ansatz computation of the susceptibility. 

5.2. Static succeptibility by the Bethe ansatz. 

Since we find a relation between lim^^o x ^ and the static susceptibility xo, we can 
perform a crucial check of our approach since %o can be computed by other means. Indeed, 
the total spin succeptibility is simply a second derivative of the free energy with respect 
to V, and the free energy can be computed directly by the Bethe-ansatz. 

By using standard manipulations for boundary theories (see eg 0]), we find for the 
impurity part of the free energy at zero temperature : 

T imp = ~ ( dd i nR (0 - 9 B )e h + {9). (5.15) 
2,11 J_ 00 av 



1 J~ ■ 

The impurity susceptibility is then given by %o = ~ t~ m ■ ^ rea ds generally 



1 / d . „, . . , 1 de%(9) 



Xo = ^>[ — lnR(A-e B )- 



g 2 7T 2 \d9 v 'V dV 



2tt d(e/2) 2 



+ /;^ lnfl(9 - M ^W), ( 5.i6) 



B=A 



where we used the correspondence e = gV. Since -4 In R(9) = cos 1 h g , we see that this is 
almost in the same form as the previous expression (|5.13| ) for xo ■ The exact correspondence 
can be established using : 



del(6) _ ( e ^(e)-Vg) de%{9) OA 
dV ~ 9+ V W dV 



V V 



from which it follows that 



de h + {9) 
dV 



29 (5.17) 
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Similarly, simple manipulations using the L operator of sections 2 and 3 lead to the iden- 
tity : 



d 2 e h + {9') yffil ^2g~ 1 dF(6'\6) 



-M8',A) 



(5.18) 

=A 



dV 2 2 V v ' ' 2 V 89 

completing the proof of the identity. 

From this discussion the static susceptibility, after a few manipulations, reads : 

xo = ( [ A de j \ P+ (e)) . (5.i9) 

ng 2 V z yj-^ cosh(6> - 9 B ) 1 J 

This expression is more convenient that the form-factors one of the previous section. In 
particular, using the Fourier representation of p, xo can be written as a convergent series. 
For small voltage this series is in V/T B (recall T B oc A 1 / 1-3 ): it corresponds to the approach 
of the IR fixed point along the stress energy tensor. For large voltage, the series is in 
(Tb/V) 1-9 , and corresponds to the (conformal) perturbation of the UV fixed point. The 
leading term in that case is found to be 

■.- 2 ( 1 -»> 1 ( T " V""" (5 20) 

By using the relation between T B and A, this can be recast as 

Xo - 2r(3 ~ 2 ^ A 2 W 2g - 3 , (5.21) 

TV 

in agreement with first order conformal perturbation theory. 

The constant T B can be easily eliminated to form universal amplitudes. For example, 
the large and low voltage results are related by the following universal product : 



hm -Xo(V)} 2(1 9) x [ lim V^x^V)] = F(g) 

Wl-g) 1 ( 2G+(i) x 2(1-9) < 5 ' 22 > 



2^(1 -g) 1 / 2G+Q 

" ^-^g^g r(-g)T(g - 1/2) cosTvg \G+(0)e A 

5.3. Singularities 

When the impurity is present, the finely tuned cancellations killing off the singularities 
will be upset by different R-matrices terms. Let us see how this works in the double well 
problem. 

Consider first uj c = gV. Recall that in the absence of impurity, there was a fine 
cancellation between the particle-hole process below uj c and the one-breather process above 
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it. In the presence of the impurity, consider first the particle-hole process. For the hole 
close to 6 = — oo, the reflexion matrix is simply equal to one. There is no dressing effect 
because the only particles in the sea affected by this hole are themselves at rapidities close 
to — oo, where their R matrix is also unity. Hence the only part that sees the impurity is 



the particle close to the Fermi rapidity. For this one we have the dressed R-matrix Q5.4T) 
so the particle hole process contributes to by a term proportional to : 



A ph { R{A) exp 



A d 



d9F(9\A) — In R(9) 
du 



1 > (u c — u)),uj < uj c . (5.23) 



At low voltage in particular this is proportional to ijr-(u) c — uj). 

The breather reflection matrix for Of, = — oo is equal to unity, so the breather process 
does not contribute linearly to x"( u ) around u c . Hence ( |5.23| ) is actually the whole leading 
behaviour close to iv c , and there is now a singularity which we refer to as being of the type 
| a; — u c \ (a discontinuity in the first order derivative). 

The same analysis can be carried out for the other processes. However, because several 
processes cross at uj c = ngV, n > 2, and there are more amplitudes than relations to fix 
them, we cannot be as conclusive: the singularity will be of the type \uj — ngV\ unless some 
special cancellations occur, in which case it will be weaker. 



6. The AC noise for tunneling in the fractional quantum Hall effect 

6.1. Low frequency behaviour 

The case of the FQHE is more complicated. In the dissipative quantum mechanics 
problem, the reflection matrices are antidiagonal in the soliton/antisoliton basis. For 
the impurity problem this is no longer true: the reflection matrix can connect a state 
consisting of all right-moving solitons to another in which there are left-moving solitons 
and antisolitons. Without a voltage, eigenstates obviously read : 

\\e > e = \e >f +R{e)i \e > £ L , . (6.1) 

At this stage, one must be very careful. We have argued before that in the presence of 
a voltage, solitons and antisolitons have different energies, so it would seem that (|6.1] ) 
is not an eigenstate of the hamiltonian when V ^ 0. This conclusion is incorrect, for a 
subtle reason. The boundary sine-Gordon model must be thought of as an anisotropic 
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Kondo model with a special (cyclic) representation of SU(2) q at the boundary |25[](or, 
more generally, a representation of the "q-oscillator algebra" |[26|| ) . The additional degrees 
of freedom provided by this boundary spin must be included in the proper definition of 
asymptotic states, and in particular they reinstate spin (charge) conservation. At the end 
of the day, this auxiliary spin can be gauged away and the result is that one can, indeed, 
treat ( |67T| ) as an eigenstate with the same energy as \9 >f . See the appendix for more 
details. 

More general eigenstates are mixtures of left and right moving particles, which we 
denote : 

||6>i, 9 n > ei ,...,e n 

In particular, the ground state gets modified, the sea of solitons becomes a sea made of su- 
perpositions of solitons and antisolitons according to \6 > + -^ \ \0 >+= \6 >+ +P(6)\6 >*1 
+Q(6)\6 > + , where we introduced the notation : 

K=(i q). (6-2) 

and the elements P, Q are given in J27j . When we compute the full current- current corre- 
lation function, different terms arise depending on the chirality of the J operators. They 
actually have different physical meanings, and we will treat them separately. 

To understand the physical meaning of the different terms, let us stress again that 
the boundary formalism is fully equivalent to a formalism with only R movers scattered 



through an impurity. The R state in (6.1) can then be considered as an "in" R-state, 



and the L states in ( |Q| ) as "out" R-states, with the same energy but possibly different 
quantum numbers. Formula ( |6.1[ ) is thus a scattering eigenstate in the traditional sense, 
and the following computations are fully equivalent to Landauer-Biittiker scattering [ESi 
as carried out in for the particular case g = 1/2. 

Let us finally emphasize that we are only concerned here with the to dependent com- 
ponent of the noise: the DC noise (which has been already computed in ||) is implicitely 
subtracted. 

The first term corresponds to JrJr. This term sees only the R-states in ( |6.1| ), and 
can be thought of as a noise purely in the imcoming channel: it is thus insensitive to the 
boundary interaction. One has, at coincident points : 

Srr = ^M- (6.3) 
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The second term corresponds to JlJr and is more complicated. It can be thought 
of as a noise between incoming and outgoing channels. Jr acting on the new ground 
state can create pairs particle-hole of solitons, create pairs solitons antisolitons, or add 
breathers. Only the first process contributes to low energy. The particle-hole pairs then 
"bounce" on the boundary where they can switch charges. They are then acted on by 
Jl and "destroyed". Observe however that in this destruction, an initial soliton can be 
replaced by an antisoliton since now the ground state is a mixture. 

The dressing follows a similar principle. If we write : 

P*(0i)P(0i) + Q*(0~i)Q(0i) « [1 + - 0i)], 



(where <p is rea l as can easily be proven using that |P| 2 + \Q\ 2 = 1) we can dress the 
reflection matrices according to : 



Q(6) = Q(0)expz 
V(6) = P(0)expi 



F{e'\e)(j){e')de' 



-oo 
A 



F(e'\e)<p(e')de' 



(6.4) 



From [|l|, we find that the low frequency noise is now written : 



S LR {u){x u x 2 ) = I dte i " t (J L {t,x 1 )J R {0,x 2 )) 

-oo 

d(3 2 d(3' 2 [V\e 2 )V{9' 2 ) - Q*(e 2 )Q(9' 2 )]e^e^x 

3 



9_ 

'2% 



(6.5) 



here X\ and x 2 physically correspond to the two sides of the impurity, and should be taken 
different in general. For the tunneling noise, both x\ and x 2 will approach and the phase 
will disappear. By the same manipulations as above this can be rewritten : 

Slr{u){x u x 2 ) = JL e M*i-*2) / [p(x)T*(x-u)-Q(x)Q*(x-uj)}dx. (6.6) 
2 ?r Jo 

Of particular interest is the limiting behaviour of S as uj — > 0. In that limit, we pick up 
the value of the argument in the integral (|6.6| ) for x = for which the dressing effect just 
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cancels out. One finds simply : 



Slr(w) « ^ [\P(A)\ 2 - \Q(A)\ 2 ] e M«i-*) 



(2|P(A)| 2 -l)e^ (xi - X2) 



2tt v ~ |J v ^ vi (6.7) 



\ e *G-VA+T B 



2(^-1) 



At large voltage it goes to the noise without impurity as expected. At small voltage we 
find : 

(,B, 

The last term corresponds to JlJl and can be interpreted as a noise purely in the 
outgoing channel. It does have quite a bit of structure. At leading order at low frequencies, 
the first process that contributes is the creation of a particle hole pair near the Fermi 
surface. Observe however that due to the mixing induced by the boundary, one can either 
destroy a soliton and create another one above the sea, or destroy an antisoliton and create 
another one above the sea 0. At leading order, one finds, at coincident points : 

S LL {u>) * ^ (\P(A)\ 2 - \Q(A)\ 2 ) 2 = 9 P(l- *\P(A)\*\Q(A)A • (6.9) 

Z7T ATT \ TV J 

Observe that Srr,Slr and Sll are respectively of order 0,2 and 4 in the R-matrix 
elements. Adding all the components we find the low frequency noise of the tunneling 
current : 

This noise reproduces the standard result St = at small voltage, where there is no 
transmitted current. This noise vanishes at very large voltage, when the tunneling current 
goes to zero, as : 



11 Due to this mixing also, it is also possible to destroy a pair of particles close to the Fermi 
surface. This process however does not contribute at leading order because the form factors for 
destroying a pair +— and a pair — h are opposite. 
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The noise for the current running down the sample can also be expressed simply as 



S T (u) « y -^\P(A)\\ (6.12) 

Let us stress here that a priori, all these results hold only for < g < 1. As g — > 1, 
the threshold for adding an antisoliton vanishes Min (e_) = (1 — g)V — > 0, so exactly at 
g = 1 other processes are implied in the low frequency physics. In other words, we expect 
a priori that 

S(u) S(iv) 
lim lim — : — j— yt: hm hm — — — . (6.13) 

Interestingly however, exactly for g = 1 the noise is easy to compute directly [p9| since 
the electrons are non interacting, and one finds the same as ( 6.11 ). In that light, ( |6.10| ) 



actually appears very natural since we showed that the low energy excitations at V 7^ 
are described by an effective free fermion theory. 

Similarly, as g — > 0, the threshold for adding breathers vanishes, and non commuta- 
tivity of the limits is also a priori expected. 

6.2. Other singularities? 

The effect of thresholds we have discussed for the double well problem should be 
observable as well in the noise, both for Slr and Sll, leading to a first singularity of the 
form \uj — gV\, and presumably other singularities of the same formm \oj — ngV\. As an 
example, let us discuss the singularity at u = gV. Like in the double well problem, it 
arises because the two processes of creation of a pair particle-hole and of a breather do not 
match at u = u c . For the particle-hole process, the bulk amplitude gets multiplied by a 
factor |Q(A)| 2 . This is because, for a hole at 9\ and a particle at 6*2, the general amplitude 
is : 

|2 



\V*(9 2 )V(e 1 )-Q*(6 2 )Q(e 



1/ 



and at leading order close to the threshold, we have 62 = —00 and 9± = A. For the 
breather process, the bulk amplitude is unchanged since at leading order, the breather 
reflection matrix appears only in the form of its modulus square. Hence the singularity is 
proportional to 

\1-Q(A)\ 2 ] \u-gV\. (6.14) 



In addition, other singularities take place in Slr. This is because, beside the processes 
of the double well problem, other processes are allowed here, for instance the destruction 
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of a pair of particles close to the Fermi surface, as already mentioned earlier. Now this 
process can take place provided the frequency ui is smaller than twice the maximum energy 
of a hole in the sea, or uj < 2gV . Hence, this process contributes a singularity which is 
however weaker than \u — 2gV\. More generally, processes where 2p particles are destroyed 
in the sea have a maximum frequency uj = 2pgV, and should lead to singularities weaker 
than \uj — 2pgV\. Recall that such processes are allowed in the general interacting theory, 
in contrast with the case g = 1/2 where only pairs can be created by the current, which 
implies that the excess noise in the outgoing channel vanishes for u > 2gV. 

6. 3. Comparison with perturbative results 



Finally, we would like to compare our results with the perturbative approach of [JTC 
1 l[ concerning the noise in the tunneling problem. For the low frequency behaviour, we 
agree that there is a singularity of the type \u\, but the amplitude ( |6.11| ) that we find does 
not agree with these authors, except for g = 1/2 and g = 1. Observe that the amplitude 
we find for the \uj\ singularity, while expanding nicely at large A in powers of A 4 / 3 , does 
not expand in powers of A 4 at small A. This suggests that the UV perturbation theory 
attempted in |T^],[^TJ has a vanishing radius of convergence. In particular, we find that 



the large voltage behaviour of the amplitude goes as l/ 4 - 4 /f, instead of l/ 4 ^" 1 ). 

In fact, we can discuss the difference in more details. Indeed, Chamon, Freed and 
Wen have recently |5(J obtained for the noise the expression 111 

2 



UJ 



dl\ 



instead of our expression fl6.12j ). It is worthwhile to examine the difference in more details. 
In the DC current in the presence of a voltage was determined, 

I = 2n [ A p + (e)\P(e)\ 2 . (6.16) 

To obtain the differential conductance, we take derivative with respect to V. The derivative 
of the density can be obtained from 



2.^ = V /§L( M ), (6.17) 



12 Unfortunately, what is called L and R here is a bit different from the notations adopted in 



[ 11 ] , [ 30 1 . The different noises are however simply related by linear combinations. 
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from which it follows that 

g = §f = \[1 \ p ( A )\ 2 + f mA)\p(e)\ 2 de 

Define a new renormalized reflection matrix 



(6.18) 



V dlff (9)\ 2 = \P(9)\ 2 - j $(9-9')\V diff (9')\ 2 d9', (6.19) 

J — oo 



then ( |6.18| ) reads 



G=\ 9 ^fM)\\ (6.20) 



The meaning of this result is as follows. Starting with a voltage V, one is interested in the 
additional current when V — ► V + 5V. This current has two origins: one is the shift of the 
Fermi sea, the other one is the change in populations deep in the Fermi sea. The shift of 
the Fermi sea adds a number of solitons 

SN = L P+ (A)5A = ^J^SV. (6.21) 

If we want to write the change of current in terms of those particles only, it has to involve 
a dressed reflection matrix which takes into account all what happens deep in the Fermi 
sea: 

^ = 2n 6 -^\V dlff (A)\^ (6.22) 

which coincides with fl6.18| ) using ( |6.19| ). Hence the differential conductance can be fully 
interpreted in terms of a new dressed reflection matrix. 
In this language, the formula (|6.15|) reads 

S T {u)^ 9 ^\V dlff {A)\ A - (6-23) 

The difference with our formula ( |6.12| ) is thus fully a difference in the dressed reflection 
matrix. Naively, based on current computations, one would have expected (eg by analogy 
with dissipative quantum mechanics) that ( |6.23| ) would hold. The key point however is 
that the dressing of reflection matrices is not a universal property, but depends on the 
quantity under study. For the conductance, the whole of the Fermi sea matters, and the 
renormalized reflection matrix V d i / / involves the whole sea as well. For the (T = 0) noise, 
the only effect of the sea in V is a phase that disappears in moduli square. Presumably, 



this is a non perturbative effect that cannot be seen in the approach of |3(J . 

Similarly, we disagree on the singularity structure. Except for g = 1/2 we find that 
the noise S t (w) has a singularity at u = gV - the "quasi particle singularity" - where the 
first derivative is discontinuous, while the authors of [ ID| , [ |11| argued that this singularity 
was either a (/-dependent power law, or was absent. 
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7. Conclusions. 



This paper is a first step towards the computation of correlation functions in quantum 
impurity problems in the presence of a voltage and a temperature. While much remains 
to be done, we think the present results already indicate very interesting features, and 
could lead to numerical and experimental applications. In particular, the presence of a 
singularity at u = gV at T = 0, while probably unobservable experimentally, should lead 
to a pic in the derivative dS }^ at finite T. We also expect that such pics should appear 
in fact at regularly spaced values u = ngV at finite T (while presumably the amplitudes 
of these pics will decrease very rapidely with n). The meaning of the singularities at 
T = can be understood in the double well problem: x"(o;) going as \u> — gV\ = \u> — e\ 
simply means that the long time behaviour of the spin spin correlator has an oscillatory 

i e t 

component -p-. 

A more complete discussion of the low energy excitations, including further justifica- 
tions of ( p.6|) , will be provided in the following paper [p2 |, together with some results at 
V = but T > 0. 
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Appendix A. The effects of the magnetic field and voltage 

Introduce : 



Lg = ttt — / dx [(d t <P) 2 - (d x <P) 2 ] . (A.l) 



Let us then define 



4n 



L(h') =L G -^— I dt<f>. (A.2) 



DO 
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Introducing 



SivgU 2 + -^(cW>) 2 
8ng 



(A.3) 



the associated hamiltonian is 



H(h') = H G + L^-(h') 2 + 2gh' [ U(x)dx. (A.4) 

More generally, we can also put a boundary coupling and a boundary field, defining 

H B = X (s+e-^ 2 + S_ e ^°V 2 ) , (A.5) 



and 



/0 7 
n(x)dx+-5 2 , (A.6) 
-oo 2 



The unconventional normalisation is taken to agree with the presentation in the bulk of 
the text. The total z-component of the spin, which commutes with the hamiltonian, reads 
with these normalizations 



S-, 



1 



total 



2ivg 



d t <j) + S z 



(A.7) 



Hence, when b! = -, ( JA.6| ) describes the field | coupled to the total spin. Otherwise, a 
different field is coupled to the boundary and bulk components of this spin. Now, the term 
linear in IT in ( |A.6|) can be eliminated by a unitary transformation. Introduce 



U = exp 



- — [ (p(x)dx 

47T 



then one finds 



9 ,ui\2 



UH(h, h^U' 1 = H(h, h' = 0)- L^-{h 



(A.8) 



(A.9) 



As a result, the thermodynamic properties of the impurity are the same whether there is a 
field coupled to the impurity spin S z only, or a field coupled to the impurity and another 
field coupled to the the component f_ dt(f>, including the particular case of a field coupled 
to the total spin. This is similar to the observations in |3l] . To be more precise, introduce 

(A.10) 



Zi(h, ti) =Tre- f,H ( h > h '\ 



with the trace taken in the spin 1/2 representation (where by convention recall we chose 
S z = ±1). Then one has 



Zx(h) 



A) 



Z x {h, h' = 0, A) 



Z 1 (h = 0, h', A = 0) Z 1 {h = 0,h' = 0, A = 0) 
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(A.11) 



In [25], we checked this relation in the particular case h! = — in the first numerator, since 



we found in |25| that the coefficients of the perturbative expansion in A of the second 
ratio in ( |A.11| ) coincide exactly with the results of TBA calculations for the Kondo model 
coupled to the total spin, ie the first ratio in ( |A.11| ). 

In QA.7 ), the g factor is a standard renormalization that can be observed for instance 



in the continuum limit of the XXZ chain regularization |2T| of the Kondo model |[32|| . For 
convenience, we define the spin of a soliton or antisoliton to be ±i in what follows. 

The scattering theory description of the spin 1/2 Kondo model involves only the 
soliton and antisoliton because it is an IR description, and in the IR the boundary spin 
is screened. Conservation of the spin then simply translates into the fact that solitons 
bounce back into antisolitons. 

A coupling \S Z to the boundary spin can also be traded for a time dependent boundary 
coupling. Indeed, consider for instance the computation of Z\ perturbatively in A. At order 
2n, one has Coulomb gas integrals related with a 2D Coulomb gas on a circle with n positive 
and n negative charges alternating. Call r = it the imaginary time, and suppose an S+ 
term has been inserted at r and an S- term at r' > r. The effect of the term %S Z in the 
hamiltonian ( |A.6| ) is to give a weight exp— (r' — r)h to the pair of insertions of S+, S-, 
relative to the weight with no insertion. Now from ( |A.1| ) we see that every insertion of 
S± is coupled with an insertion of exp=pz(/>(0)/2. The effect of the ^S z term can thus be 
absorbed into a time-dependent phase multiplying the vertex operators in Hb, 

S± -> S±e ±uh . (A.12) 

In addition, the term \S Z gives rise to an overall term exp±/3/i/2, independent of the 
order n, that can be absorbed by taking a modified trace Tr(.) — > Tr(.) e^ Sz . 

The foregoing chain of arguments generalizes to higher spins in SU(2) q . This is 
obvious for the h! term which is independent of the spin. As for the h term, it follows 
simply from the fact that insertions of S± change the spin S z by ±1, independently of the 
representation. 

The properties of correlators with ( |A.6| ) can be analyzed by the same sort of arguments, 
and the various transformations easily followed on each term of the correlators. 

The first relevant point for this paper is that, although the dissipative quantum me- 
chanics (anisotropic Kondo) hamiltonian has a field coupled only to the boundary spin 
h = e = gV, it can be transformed, up to a shift of the ground state energy, into a hamil- 
tonian where the field is coupled to the total spin h' = - = V, h = e = gV, and then easily 
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diagonalized by Bethe-ansatz. In the massless description we are using, the boundary spin 
is screened, and only solitons and antisolitons have to be considered, with an energy shift 
equal to ±V/2 = ±e/2g. Since we are using a scattering approach, the eigenstates are built 
as explained at the beginning of section 4 by combining the left and right (asymptotic) 
eigenstates of the full line hamiltonian ( |2.1|) . The presence of the impurity does actually 
change the densities by a factor of order 1/L, which does not affect the results we are 
interested in. 

The second point concerns the impurity in a Luttinger liquid. Recall, following ||25|| , 
that the boundary sine-Gordon model is in fact equivalent to an anisotropic Kondo model 
where the boundary spin is defined in a special, cyclic representation of SU{2) q (or, alter- 
natively, in a representation of the q-oscillator algebra |33|). While these boundary spins 
can be "gauged away", let us discuss them a bit more. First, asymptotic states are now 
characterized by particles and by the value of the renormalized boundary spin. Indeed, 
let us stress here that the boundary spins that appear in the hamiltonian sh and in the 
scattering theory s sca t differ for two reasons. There is the screening, already mentioned 
previously for the Kondo problem, and the multiplicative renormalization, which can be 
observed by considering the energies with a magnetic field. One finds s sca t = g( s H — 1), 
where the g factor was also mentioned previously. 

The main use of the boundary spins is to reinstall conservation of the spin. This is 
obvious on the hamiltonian which looks then like a higher spin Kondo problem. This is also 
obvious in the scattering theory : every time a soliton bounces back as a soliton instead 
of an antisoliton, the boundary spin increases its value accordingly. The conservation of 
the spin allows diagonalization of the problem with a voltage. In particular since a pair 

R-soliton, S z = — and L-soliton, S z = have the same charge and the same kinetic 

' z g ' z g to 

energy, they should have the same energy, ie e_ — ^ = e + + "j 5 where we used h = gV , an 
identity which we mentioned previously. Also, in section 4, while the state (|6.1|) seemed 
to be made of particles reflected with different energies, it is more precisely a combination 
of true asymptotic states that involve also the boundary spin, that all have the same total 
spin, and all the same energy. Therefore, for the impurity problem, one starts with a time 
dependent term in the boundary coupling in ( |1 . 2| ) . This term is then transformed into a 
field h = gV applied to the boundary spin. An additional term coupled to the component 
f_ dt<p is then added, and allows for the use of the Bethe ansatz with a field coupled to 
the total spin. The boundary spin is finally gauged away. 
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The arguments of this appendix provide the missing steps in the proof of the DC noise 
formula (27) of ||. They also justify why in [25| the Bethe ansatz for a field coupled to the 
total spin was used, although the initial hamiltonian had a field coupled to the boundary 
spin only. 

Finally, let us stress that, while the double well is a problem in equilibrium, the 
tunneling problem is not. As simple way to see that is to observe that, while we used an 
eigenstate of the boundary hamiltonian with the R part made only of solitons, there are, 
thanks to the boundary spin which can absorb all spin (charge) shifts, an infinity of other 
boundary states, and averaging over all of them would, for instance, lead to a vanishing 
DC current as expected in equilibrium. Our approach here is fact fully equivalent to the 
Landauer-Biittiker scattering. In that context, it is also interesting to mention that the 
transformations previously discussed change the boundary conditions at infinity, which 
can be interpreted in terms of reservoirs. The latter are implicitly there in the scattering 
approach [|TTJ . 
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